Ground state of three qubits coupled to a harmonic oscillator with ultrastrong 

coupling 
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We study the Rabi model composed of three qubits coupled to a harmonic oscillator without 
involving the rotating-wave approximation. We show that the ground state of the three-qubit 
Rabi model can be analytically treated by using the transformation method, and the transformed 
ground state agrees well with the exactly numerical simulation under a wide range of qubit-oscillator 
coupling strengths for different detunings. We use the pairwise entanglement to characterize the 
ground state entanglement between any two qubits and show that it has an approximately quadratic 
dependence on the qubit-oscillator coupling strength. Interestingly, we find that there is no qubit- 
qubit entanglement for the ground state if the qubit-oscillator coupling strength is large enough. 

PACS numbers: 42.50.Ct, 42.50.Pq, 03.65.Ud 



I. INTRODUCTION 

Recent experimental progresses related to the qubit- 
oscillator system in ultrastrong coupling regime have 
been re por ted in different light-matter interaction sys- 
tems [iHlOl , where the coupling strength between a single 
qubit and a single oscillator reaches a significant fraction 
of the oscillator and qubit frequencies. In this ultrastrong 
coupling regime, the ubiquitous Jaynes-Cummings model 
[Tlj under the rotating-wave approximation (RWA) is ex- 
pected to break down leading to a mass of unexplored 
physics and giving rise to unconventional quantum phe- 
nomena [12-2^]. For example, superradiance transition 
[13] , vacuum Rabi-splitting [l], , photon blockade [13] , 
Bloch-Siegert shift and plasmonic effect [Toj . 

Since the Hamiltonian of a qubit-oscillator system con- 
tains counter-rotating terms that make the computa- 
tional subspace unclosed, the fully analytical solution to 
the ground state of this Hamiltonian in the ultrastrong 
coupling limit is still not found. Although the spectrum 
and eigenfunction of the Rabi model beyond the RWA are 
known by numerical diagonalization in a truncated finite- 
dimensional Hilbert space [H,[lll, the analytical solution 
to the qubit-oscillator system beyond the RWA is quite 
necessary for capturing the fundamental physics clearly. 
Such an analytical treatment has the potential to be ex- 
tended to more complicated models for the implementa- 
tion of the quantum information processing (QIP) [27j . 
Therefore, various mathematics approaches have been 
proposed to analytically obtain the ground state prop- 
erties of the sing le-qubi t Rabi model in the ultrastrong 
coupling regime 



RWA mehtod ^0 



|. For example, the generalized- 
|3l| functions well when the qubit fre- 
quency is smaller than the oscillator frequency, the varia- 
tional treatment for the ground state ^,,34| captures 
reasonably in the single-qubit Rabi system which is very 
hard to be generalized to multi-qubit Rabi systems, and 
the transformation method is a perturbation expansion 
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and has been successfully applied to the single-qubit Rabi 
system [33,li3-Iil- 

Recently, the Tavis-Cummings model beyond the RWA 
has been extended to the multi-qubit case using an adi- 
abatic approximation method when the qubit frequency 
is far larger than the oscillator frequency ^1], and the 
ground state of the nearly resonant Rabi model of two 
qubits coupled to a harmonic oscillator has been ana- 
lytically treated by usin g b oth the variational and the 
transformation methods [4J|. The Rabi model of three 
and more qubits coupled to a common oscillator in the 
ultrastrong coupling regime has more potential applica- 
tions in QIP [2l|, [23] than the single-qubit Rabi model, 
such as protected quantum computation [21], which is 
expected to be very promising with the circuit QED ar- 
chitecture. 

However, the ground states of the three- and more- 
qubit Rabi models in the ultrastrong coupling regime 
have not been extensively studied. Braak generalized the 
method based on the Z2 symmetry [s^ to the three-qubit 
Dicke model [13] to analytically determine the system's 
spectrum, which is dependent on the composite transcen- 
dental function defined through its power series. How- 
ever, this method can not be extended to determine the 
concrete form of the ground state. 

Different from the Ref. [l^l , we focus here on the ana- 
lytic ground state of the three-qubit Rabi model in the ul- 
trastrong coupling regime by the transformation method. 
By mapping the three-qubit Rabi model into a solvable 
Jaynes-Cummings-like model, we show that the ground 
state energy and the ground state of this three-qubit Rabi 
model can be approximately determined by the analytic 
expression based on the transformation method, which 
agrees well with the exactly numerical simulation in the 
ultrastrong coupling regime under different detunings. 
The ground state entanglement between any two qubits 
is characterized by using the pairwise entanglement and 
has a quadratic dependence on the qubit-oscillator cou- 
pling strength, which is able to be analytically deter- 
mined within a wide range of parameters. The inter- 
esting feature in the ground state entanglement exists 
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in its maximum value, which decreases quickly to zero 
and never increases again as the qubit-oscillator coupling 
strength is large enough. 



II. TRANSFORMED GROUND STATE 

The Hamiltonian of three identical qubits coupled to 
a harmonic oscillator without the rotating-wave approx- 
imation is {h — 1) 



H 



1 



Wa{J+ + J-) + Wca^a + g{a'' + a)Jz, (1) 



where a and are respectively the annihilation and cre- 
ation operators of the harmonic oscillator with frequency 
Wc- Ji {I ^ ±,z} describes the collective atomic op- 
erator of a spin- 1 system, satisfying the angular mo- 
mentum commutation relations [Jz,J±] — ±t/± and 
[J+, J_] = 2Jz. Physically, the spin-| system is non- 
trivial and the states are entangled in terms of individ- 
ual qubit configurations. Wa is the transition frequency 
of each qubit. g represents the collective qubit-oscillator 
coupling strength. 

The key point in this paper is to determine the ground 
state energy Eg and the ground state \(j)g) for the three- 
qubit Rabi system in the ultrastrong coupling regime, 
where H\(t)g) = Eg\(t)g). To derive the analytic ground 
state, we define \m)a to be an eigenvector of Jz, i.e., 
Jz\m)a = m\m)a {m = — §,— ^,^,1). Besides, we will 
respectively use |^)/ and |0)/ to represent the coherent 
field state with the real amplitude X and the vacuum 
field state. 

In what follows we extend the transformation method 
used in the single- and two-qubit Rabi models ^] to 
the three-qubit Rabi model. To transform the Hamilto- 
nian H into a mathematical form without the counter- 
rotating wave terms, we apply a unitary transformation 
to the Hamiltonian H: 



with 



(2) 



(3) 



where x is a variable to be determined. Then the trans- 
formed Hamiltonian H is decomposed into three parts 
[431: 



H — Hq + + H2 



(4) 



with 



i?o = iTu^aJx - (2.gx - WcX^)JI + WcO^a, (5) 

H^ = {g- WcX){o) + a)Jz + irjWaXia^ " a)Jv, (6) 
H2 = WaJx{cosh [x(a^ - a)] - r/} 

+iWaJy ^sinh [x{a^ - a)] - 77x(a^ - a)|, (7) 



where rj =/(0| cosh[x(at - a)]|0)/ = exp[-xV2]. The 
terms cosh [x(a^ — a)] and sinh [x(a^ — a)] in H2 have 
the dominating expansions: 



cosh[x(at - a)] = r] + 0{x^), 
sinh[x(a''' - a)] = xvia'^ - a) + O(x^), 



(8) 
(9) 

here O(x^) and 0{x^) represent the double- and multi- 
photon transition processes containing higher-order op- 
erators for a^ and a, which can be neglected as an ap- 
proximation when X is small. Thus, H ~ Hq + H^. By 
now, ou r ap proximation procedure is the same as that 
in Ref. [43]. However, the main difference exists in the 
diagonalization for with the collective qubit vectors. 
Such an operator Hj — rjWaJx — (2(?x ~ WcX^)J'^ appear- 
ing in the qubit basis Fa = {| - |)a, | - 5)0, |^)a, ||)a} 
represents a renormalized four-level qubit system. This 
is different from the two-qubit Rabi system [43| which is 
a renormalized three-level atomic system corresponding 
to the diagonalization for Hq. Through diagonalizing the 
operator Hj in the qubit basis Ta, we obtain the renor- 
malized qubit eigenvector \ipk)a with the eigenvalue Afc 
{k = 1, 2, 3, 4) as following: 
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and 
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where iV^ = \/2 + 2K'^ (k = 1,2,3,4) is the normaliza- 
tion factor for the eigenvector \ipk)a- For x ~ 5 ~ 
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FIG. 1. (Color online) (a) Numerical solutions for the variable 
X that makes Ci = 0. (b) The C3 value as a function of g 
for different qubit-oscillator detunings. (c) The ground state 
energy as a function of the coupling strength g. The solid 
line represents the transformed ground state energy Eg and 
the dash line represents the exact ground state energy Eg. 
(d) The fidelity F of the ground state \(j)g) obtained by the 
transformation method. 



-T]WaX [ViK^ + 2K1K2 - V3K2^ , (13) 

C3 = {3 + KiKi)ig - WcX) 

-VWaX (^VSKi + 2K1K4 - VSKi^ . (14) 

Therefore, |(/9i)a|0)/ is expected to be the approximate 
ground state vector if the conditions Ci = and C3 ~ 
are both satisfied, then the ground state \(j)g) of this 
three-qubit Rabi system approximates the transformed 
ground state \<j)g)- 



11 3 3 

-Kl\-)a\ - -X)f + \^)a\ - -X)f 



(15) 



and the transformed ground state energy E is: 



the eigenvalues here are arranged in the decreasing or- 
der through the numerical simulation, i.e., Ai < A2 < 
A3 < A4. Therefore, H can be expanded with the above 
renormahzed eigenvectors: 



fc=l 



+ {C3a + C4a'')\ipi)a{^4\ + (Csa + C6a^)\(p2)a{(P3\ 



+ iCra + Csa^)\ip3)a{v4 



H.c. 



+ WcCi a, 



(12) 



where Cx{x — 1, 2, 3, 8) is the coefficient depending on 
the variable x- It is obvious to see that Ci,C3,C5 and 
C7 represent the coupling strengths of the corresponding 
counter-rotating terms with respect to the renormahzed 
eigenvectors in Eq. (|12p . 

Similar to the single- and two-qubit Rabi systems 
(4^ I43I , the main task after transforming the Hamilto- 
nian H into H is to eliminate the counter-rotating terms 
for the eigenvector with the lowest eigenenergy. The ma- 
jor obstacle here is to remove the different coupling coef- 
ficients Ci and C3 of two counter-rotating terms for the 
eigenvector \ipi)a simultaneously. This is very different 
from the single-qubit .42.] and the two-qubit Rabi mod- 
els [13 , which just have one counter-rotating term for the 
approximate ground state vector. Although it is not pos- 
sible to simultaneously remove the coefficients Ci and C3 
for all the values of x, we find that the conditions Ci — 
and C3 « can be both satisfied when < x ^ 0.5, 
meaning two counter-rotating terms for the approximate 
ground state vector \ipi)a can both be eliminated if the 
qubit-oscillator interaction is not too strong. The coeffi- 
cients Ci and C3 have the following analytical forms: 

Ci = (3 + KiK2)ig - w^x) 
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According to the condition Ci = 0, the numerical solu- 
tion of X is plotted as a function of the coupling strength 
(7 for different qubit-oscillator detunings in Fig. 1(a). We 
find X has a proportional relation with g: x — w 4-^, ■ 
By substituting the result X from Fig. 1(a) into Eq. (14), 
we obtain the corresponding solution for C3 in Fig. 1(b), 
which shows the conditions Ci = and C3 ~ can be 
both satisfied when < g < Q.'awa- This guarantees 
two counter- rotating terms in the eigenvector \^pi)a are 
both eliminated if the qubit-oscillator coupling is not too 
strong. In Fig. 1(c), we have estimated the accuracy 
between the transformed ground state energy and 
the exact ground state energy Eg under different qubit- 
oscillator detunings, in which Eg has an approximately 
quadratic dependence on g: 



E„ 



2Wc + iWa 
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(17) 



We see that the transformed ground state energy achieves 
the nearly perfect matching with the exactly numerical 
value within the ultrastrong coupling regime g < 0.5wa- 
When g = O.bwa, the errors for the transformed ground 
state energy at Wc = O.Swa, Wc = Wa, and Wc = 1.2wa 
are 0.49%, 0.19%, and 0.07%, respectively. Especially, 
when there is a positive detuning Wc — Wa > 0, the trans- 
formed ground state energy fits much better with the ex- 
act value for a wide range of g than that with the negative 
qubit-oscillator detuning or the exact qubit-oscillator res- 
onance, and its error is only 0.9% even when g = O.Swq 
for Wc = 1.2wa- This result coincides with the variational 
behavior of C3 in Fig. 1(b), in which C3 grows much 
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FIG. 2. (Color online) The pairwise entanglement Np^ as 
a function of the coupling strength g under different qubit- 
oscillator detunings: (a) Wc — O.Swa; (b) Wc = Wa; (c) Wc = 
1.2wa. The solid line (dashed line) corresponds to the pairwise 
entanglement of the transformed (exact) ground state. 



slower with Wc = 1.2wo than the case with Wc = O.Swa or 
Wc = Wa when the coupling strength satisfies g > O.Swq. 

To examine the reliability of the transformed ground 
state \4>g'), we use the fidelity F, which is defined as 
F = {(pg\(j)g), as a measurement for the transformation 
method. From the result of Fig. 1(d), we find that the 
exact ground state for the three-qubit Rabi model can be 
approximately represented by \(j)g) within the ultrastrong 
coupling regime < g < 0.5wa- For example, we obtain 
the fidelity with high value F > 99% for g < O.bWa under 
different qubit-oscillator detunings. 



III. GROUND STATE ENTANGLEMENT 

To investigate the qubit-qubit entanglement for the 
present three-qubit Rabi model in the ground state, in 
which the prescription set out for symmetric Dicke states 
is used, we proceed to consider the pairwise entanglement 
[45, 46] between any two qubits. 

Taking the transformed ground state \(j)g) in Eq. ([T5|. 
the reduced density matrix pa of any two qubits can be 
written as: 



Pa 



^ pii pi4 

P22 P23 

I P32 P33 

\ P4,i pm 



(18) 



where 

Pii = Paa 

Pii = Pii 



-2N + 4( J2) 1 



miN - 1) 6 3(1 + K^) ' 



N{N-l) i{l + Kl) 

- 4( J2) _ Kl 



P22 — P23 — P32 — P33 — 



AN{N~l) 3{1 + Kf)'' 



(19) 



and the standard basis in pa is { \ei)\em), \ei)\sm)^ 
) }, with \ei) (|e™)) and \si) (|sm)) 



{l,m = 1,2,3; and I ^ m) denoting the excited and 
ground state of the Ith (mth) qubit, respectively. There- 
fore, the pairwise entanglement Np^ can be expressed as: 



N„ 



2 max <^ 0, \p23\ - y/PllPM, \pu\ - y/P22P33 



(20) 



In the ultrastrong coupling regime g < 0.5wa, we can 
numerically verify: 



Kl) 



4(M;a + Wc) 



r5'- (21) 



Fig. 2 illustrates the pairwise entanglement Np^ obtained 
from the transformed and the exact ground states versus 
the coupling strength g under different detunings. We 
see that the pairwise entanglement has a quadratic de- 
pendence on g at small coupling strength, which is math- 
ematically captured by the approximate power law be- 
tween Np^ and g in Eq. (pOj) . If g > 0.5wa, discrepan- 
cies for the numerical results between the transformed 
and exact ground states become bigger as the coupling 
strength increases further. The maximal pairwise entan- 
glement between any two qubits is determined by the 
detuning IS. — Wc — Wa- For the positive detuning A > 0, 
the increase of A leads to the increase of the maximal en- 
tanglement in the qubit system, whereas it is oppositive 
for the negative detuning case. Interestingly, the pair- 
wise entanglement Np^ decreases quickly to zero after 
reaching its maximum, and remains at zero even when 
the coupling strength g increases further, which means 
that there is no qubit-qubit entanglement in the ground 
state of such a model any more if the coupling strength g 
is large enough. For example, the pairwise entanglement 
decreases to zero for g = l.bwa at the exact resonant case 
and never increases again even when g increases further. 
This feature is distinguished from the result of Ref. [11] 
and can be explained as follows. When the field and one 
qubit is traced out, the first and fourth terms of Eq. (15) 
do not result in the entanglement of the other two qubits. 
In other words, the pairwise entanglement is contributed 
by the W-si&te components \ — \)a and W)a- The coeffi- 
cient Kl of the terms involving \ — \)a and \\)a quickly 
decreases to zero when g is large enough resulting in the 
vanish of pairwise entanglement. 



IV. CONCLUSION 

In summary, we have shown that the ground state 
of the three-qubit Rabi model in the ultrastrong cou- 
pling regime can be analytically treated by the trans- 
formation method. The transformed ground state fits 
very well with the exact ground state for different de- 
tunings even when the coupling strength g increases to 
0.5z«a. When g = 0.5u'a, the error of the transformed 
ground state energy is only 0.19% at Wa = w^, and the fi- 
delity for the transformed ground state keeps higher than 
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99% when g < 0.5wa under different qubit-oscillator de- 
tunings. Finally, we use the pairwise entanglement to 
analytically examine the qubit-qubit entanglement, and 
the result shows that the ground state entanglement has 
an approximately quadratic dependence on the qubit- 
oscillator coupling. Interestingly, we find that there is 
no ground state entanglement if the qubit-oscillator cou- 
pling strength is large enough. 
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